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Abstract 

We recalculate the hadronic contribution to the effective electromagnetic 
coupling constant a(M^) using quasi-analytical approximations for the cross- 
section of the e~^e~ — > hadrons proposed earlier by one of the authors (N.V.K., 
Mod.Phys.Lett.A9(1994)2825). We find that a-i(M|) = 128.98 ± 0.06. 
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For LEPl observables at the Z-pole the largest effect from pure QED corrections 
is the change in the effective electromagnetic constant when going from = 0, 
where the fine structure constant a^^ = 137.036 is measured to = M| The 
change from a = a{q^ = 0) to a(M|) is related to the photon vacuum polarisation 
function n(g^) via the relation 

where in the leading log approximation 

n(g2) = n;epi(g') + n;,,,(g2), (2) 

n,..tf)^iii:(>n(^)-^0(^)) (3) 

and the hadronic contribution UhadiQ"^) can be directly determined by the total 
cross-section of e^e~-annihilation into hadrons [|1| 

Re{UUs)) = '^P r , ds\ (4) 

37r J4m2 s [s — S) 

ah{s) = crtot(e+e" ^ 7 hadrons) = — — R{s) (5) 

oS 

So the main problem in the calculation of a(M|) is the evaluation of the integral 
(4). There are several papers H] - P devoted to the calculation of the integral (4). 

In this paper we recalculate the hadronic contribution to a(M|) using quasiana- 
lytical approximations for the cross-section of e^e^ hadrons proposed earlier by 
one of the authors |Q. The motivation for such recalculation is the appearance of 
new data on e+e~ - annihilation into hadrons in low energy region and new value of 
the effective strong coupling constant a;s(M|). We find that 

a-i(M|) = 128.98 ±0.06 (6) 

In ref.[^ the integral (4) has been calculated using quasianalytical approximation 
for R{s). Namely, low energy contribution for R{s) has been taken from experi- 
mental data and the theoretical ansatz has been used for high energy contribution 



for R{s). Two slightly different methods have been used. The first method nu- 
merically accounts for all resonances {p, u, 0, J/ip, •••) and low energy continuum 
region ^/s < 2.3 GeV. The continuum contribution to R{s) is determined by step- 
like approximation. In the approximation when strong coupling constant = 
our ansatz has the form [^] 

R%s) = 26(3 - si) + ^e{s - S2) + ^9{s - S3), (7) 

where we took si = {2.3GeVY, S2 = {4:AGeVy and S3 = {12GeVf. The parameters 
S2 and S3 correspond to the charm and beauty thresholds correspondingly. So in the 
first method we use the relation 

Uhad{s) = li2m^-2.3Gev{s) + Urj/^^r{s) + Ucis) (8) 

Here ^2m^-2.3Gev{s) is the contribution of the low energy region a/s < 2.3GeV, 
n^^j/^ x is the contribution of J/ip- , T-resonances and their radial excitations and 
Ilc{s) is the contribution of the continuum which is described by the formula (7) in 
zero approximation. In our paper we use the value 

^2„^.-2.3Gey(M|) = (6.06 ± 0.25) • 10-=^ (9) 

for the contribution of the low energy region. For the J/ip-, T-resonances and their 
radial excitations we use the formula 

n...(M|)^i:^(^. (10) 

where Mj and T^e^i are the mass and leptonic width of the i-th resonance, respectively, 
and the effective QED coupling constant at the resonance scale is used. An account 
of QCD and quark mass corrections lead to the appearance of the factors in formula 
(7) 
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The coefficients /2(0) and /3(0) have been calculated in refs. |T0[ and the function 
fi{-^) is approximately determined by the expression 



v{x) = ^{l-Ax) (13) 

The second method of the calculations consists of taking into account low energy 
region -^i < 2.3 GeV in the calculation of the integral (4) and an account of the 
high energy contribution to Tihadi^z) is performed by taking into account nonzero 
c- and b-quark masses. So in the second method our ansatz has the form 

Ilhad{s) = Il2m^-2.3Gevis) + U,{s), (14) 

where in the approximation when strong coupling constant = our ansatz for 
Rcontis) is 



6 S ^ S 6 S \ S 

(15) 

Nonzero c- and b-quark masses in formula (15) effectively take into account the 
contribution of J/ip- and T-resonances and their radial excitations. Such approach 
works rather well in the method of QCD sum rules. 

In our numerical calculations we used the value of the effective strong coupling 



constant equal to |11 



a,(M|) = 0.119 ±0.006 (16) 



and the values of the c- and the b- pole quark masses equal to |T2 

^poie ^ (1 ^ 4 a,(me(m,)) ^ ^ ^ ^ . ^ 0.1)GeV, (17) 

3 71 

= (1 + 1^£(!!!^(!M + ...)m,{m,) = 4.5 • (1 ± Q.l)GeV (18) 

3 TT 



In the method 1 the inverse effective electromagnetic couphng constant a ^(M|) 
is represented in the form 

a-i(M|) = a-^-Aa-^(0-Aa-i(v^ < 22GeV)-^a-^{J j^, T)-^or^{c)-^or^(t) 

(19) 

Here Aa;^^(/) is the leptons contribution, Aa~^(y^ < 2.3G'e\^) is the low energy 
contribution, ^a~^{J/il), T) is the contribution of J/ip--, T-resonances and their ra- 
dial excitations, Aa;^^(c) is the continuum contribution and Aa^^(t) is the top-quark 
contribution. In the estimation of Aq;^^(/) we have used two loop approximation. 
To estimate the uncertainties we assumed the uncertainties related with the choise 
of S2 and S3 to be equal 20 pecent. In our calculations we used the values f2{-j-) 
and fsi^) at = 0. However the dependence of our results in method 1 on 
the value of c- and b-quark masses is rather small and also three and four loop 
contributions are also small numerically. As it has been mentioned before in the 
estimation of low energy hadron contribution into a~^(M|) we used the results of 
ref . . Numerically we have found [| 

A«"^(/) =4.313, (20) 

Aa-^(v^< 2.3Ge1/) = 0.830 ±0.034, (21) 

Aa-\j/ij, T) = 0.160 ± 0.008(0.016), (22) 

Aa-\c) = 2.757 ± 0.019(a,) ± 0.024(/i.c.) ± 0.003(mc, m^) ± 0.032(^2) ± 0.008(^3), 

(23) 

4 1 

Aa~\mt = 17AGeV) = f -0.008 (24) 

3n 15 rrit 

In formula (23) the uncertainties 0.019(as), 0.003(mc, m^), 0.032(s2), 0. 008(53) are 
determined by the uncertainties in the determination of as(M|), c- and b-quark 
^In the calculation of Aa~^{J/^, T) we have used the data from ref. Q 
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masses, parameters of the spectrum S2 and S3. The uncertainty 0.024(/i.c.) is our 
estimate of the higher order contributions. We assume that such uncertainty is equal 
to one halph of the difference between calculated value which takes into account 
QCD corrections up to 3 loops and one loop contribution or numerically it coincides 
with 3 loop contribution. Note that we assumed 20 percent uncertainty in the 
determination of S2 and S3 that is rather conservative estimate. In formula (22) in 
the estimation of the overall error we assumed the errors in the determination of 
the contribution of different resonances are statistically independent and obtained 
the error equal to 0.008. The number in brackets corresponds to the case when we 
simply sum up the errors from different resonances. Assuming that all errors are 
statistically independent we find 



In formula (25) the number in brackets corresponds to the case when we simply sum 
up the errors. 

In method 2 the inverse effective electromagnetic coupling constant q;~^(M|) can 
be represented in the form 

a-\Ml) = a-^ - Aa-\l) - Aa-\^/^ < 2.3GeV) - Aa'^c) - Aa'^t) (26) 

Here Aa~^{c) is the contribution of heavy quark continuum and heavy quark reso- 
nances. Numerically we have found 

Aa-\c) = 2.950 ± 0.033(^3) ± 0.038(/i.c.) ± 0.030(mc) ± 0.009(mb) (27) 

Here uncertainties 0.033(0^), 0.038(/i.c.), 0.0030(mc), 0.008(mb) are the uncertain- 
ties determined by the uncertainties of Q!s(M|), higher order corrections, c-quark 
mass an b-quark mass. Assuming that all errors are statistically independent we 




(25) 



find that 




(28) 
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in the method 2. The number in brackets in the formula (28) corresponds to the 
simple summation of the errors. So we have found that both methods 1 and 2 
lead to similar values for a;~^(M|) with the similar errors. However we believe that 
method 1 is more reliable since the errors related with the as(M|) uncertainty and 
the uncertainty of higher order corrections are smaller in method 1. Besides, in the 
estimation of the errors in method 1 we assumed very conservative estimates in the 
uncertainties related with the choise of S2 and S3 (20 percent). For instance, if we 
assume 10 percent uncertainty in the choise of S2 and S3 our error in formula (28) 
will be 0.05(0.10). Therefore we quote in the abstract our estimate of a~^(M|) 
obtained using the method 1. The value obtained in our paper is very similar to the 
value a"HM|) = 128.97±0.06(exp.) ±0.07(t/ieor.) obtained in ref.@]. The decrease 
of the errors is related mainly with the better calculation of low energy contribution 
Aa~^{M'^){^/s < 2.3GeV) and the better determination of the as{M^). As it has 
been mentioned before there are several recent calculations of a~^(M|). Table 1 
shows a comparison of some recent estimates of a~^(M|). 

To conclude, in this paper we have recalculated the value of a^^(M|) using two 
methods of ref.[0 and new experimental data. We have found that both methods 
give similar results, however as it has been mentioned before we believe that the 
first method is more reliable. 

We are indebted to RFFI-DFG research program project No. 436 RUS 113/227/0 
which made possible our collaboration. 
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Table 1. Some recent values of a ^(M|) 



a-\Ml) 


ref. 


128.87 ±0.12 


[3] 


128.97 ± 0.06(exp.) ± 0.07(i/ieor.) 


[4] 


128.99 ±0.06 


[5] 


128.89 ±0.06 


[6] 


128.96 ±0.06 


[7] 


128.896 ± 0.090 


[8] 


128.98 ±0.06 


this paper 
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